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ABSTRACT 


The ice impact problem on the engine blade made of layered composite is simulated. The 
ice piece is modeled as an equivalent spherical object and has the velocity opposite to that o 
the aircraft with direction parallel to the engine axis. Near the impact region and along the 
leading edge, the blade is assumed to be fully stressed and undergoes large deflection. A 
specified portion of the blade around the impact region is modeled. The effect of ice size and 
velocity on the average leading edge strain are investigated for a modified SR-2 model 
unswept composite propfan blade. Parametric studies are performed to study the response 
due to ice impact at various locations along the span. Also, the effects of engine speed on 
the strain and impact displacements are discussed. It is found that for a given engine speed, 
a critical ice speed exists that corresponds to the maximum strain and this critical speed 
increases with increase in the engine speed. 


INTRODUCTION 


At high altitudes, when aircraft flies through clouds of super-cooled water droplets, ice 
formation occurs on forward facing structural components. One such component is the 
engine inlet. With time, the ice accretes on the inlet and eventually sheds due to structural 
vibrations. A schematic of this phenomenon is shown in Figure 1. As a result, blocks of ice 
travelling at the speed of aircraft impacts the engine blades rotating at high RPM. This 
process may cause severe damage to the blade and subsequently to the engine. In order for 
the blade to sustain the ice impact, it is necessary to properly account for these constraints 
during the design. 


Fibrous composites are ideal for structural applications such as high performance aircraft 
engine blades where high strength-to-weight and stiffness-to-weight ratios are required. 
These factors along with the flexibility to select the composite layup and to favorably orient 
fiber directions help limit the impact damage and stresses arising from large rotational 


speeds. 
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The objective of this paper is to simulate ice impact on an engine blade. SR-2 unswept 
composite propfan blade is considered for local damage analysis. The impact analysis is 
carried out by modifying the foreign object damage option provided in Reference [1]. It is 
assumed that the damage is severe in the local region and hence only a specified porticn of 
the blade around the impact region is modeled. Furthermore, large scale deflections 
accompany the impact event in the impact zone [2]. In order to simulate this behavior near 
the impact region and along the leading edge, the blade is assumed to be fully stressed e nd 
undergoes large deflection. This is accomplished in the code by modifying the membrane 
stiffness in the radial direction to reflect a fully yielded condition. A detailed justification 
for these assumptions is also given in Reference [2], Parametric studies are performed lo 
study the response of the blade due to ice impact at various locations along the span. The 
effect of ice size and velocity on the average leading edge strain and maximum impact 
displacements are investigated. 


IMPACT ANALYSIS 


Geometry oj_ Ice Impact 

The geometry of ice impact on the leading edge of the blade is shown in Figures 2(a,b). 
Impact velocity direction relative to the blade is a function of aircraft speed and the 
rotational speed of the blade (Figure 2a). The resulting impact force not only depends on 
the magnitude and direction of relative velocity but also on the mass of ice piece. 
Depending on the spacing (i.e., number) of blades and relative velocity as shown in Figure 
2b, only a portion of the ice piece hits the blade . The ice piece that is approaching the blade 
under consideration is sheared off by the adjacent blade and only a part of it impacts the 
leading edge. Effectively, the size of the ice piece that finally impacts the blade depends on 
the blade spacing, blade speed and aircraft speed. 

Modeling of the Blade 

The damage due to impact is considered to be highly localized and hence only a portion of 
the blade around the impact region (i. e., a specified portion along the span and half of the 
blade along the chord as shown in Figure 3) is modeled. In the code, the blade geometry is 
input in the form of finite element grid and nodal thicknesses. The spanwise impact region 
is specified with two parameters; namely, lower and upper bounds of radial fractions, a and 
b. However, the modeled impact span is approximated to be the region between the two 
finite element radial stations that are defined in the blade geometry and are nearest to a and 
b. Along the chord, only half of the blade is considered for impact analysis. The finite 
element used is similar to the NASTRAN (TRIA3) three node triangular plate element 1 3,4]. 
A total of 35 nodes and 48 elements are used (Figure 3). The impact is considered at the 
midpoint of the local patch (finite element node 16) along the leading edge. All the edges 
except the leading edge are considered fixed as the impact response beyond these boundaries 
is negligible. 
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Large scale deflections usually accompany the impact event. This phenomenon is simulate 
in the elements that are close to the impact node by modifying the membrane stiffness in 
the radial direction (to reflect fully yielded condition) and zeroing-out the spanwise bending 
stiffness. These special elements are shown as shaded in the Figure 3. If O yield is the 
effective yield stress of the element in the radial direction; w 2 and w 3 are out of plane 
displacements at 3 nodes; the modified stiffness associated with these degrees of freedom is 

given by [2], 


[K] — ~r Oyield te 

2x 2 


1 -1 0 

-1 1 0 

0 0 0 


( 1 ) 


where t e is the element thickness and x 2 and y 3 are the element local coordinates (Figure 4). 
The terms of this stiffness matrix replace the spanwise bending stiffnesses of the linear 
elastic triangle. The stiffness in Equation (1) reflects the perfectly plastic condition with 

constant stress (CT yie id) in the radial direction. 


Figure 5 shows the material configuration of the composite blade. In each finite element, 
the skin (material 2) is always present. The remaining materials exist only if the blade 
thickness permits them. Materials (n-l), in-2 ) ,.. are dropped from the layup of an element 
in that order. If the element thickness can accommodate all the materials then the core (n 
material with variable thickness) is used to fill the remaining excess thickness. 


Modeling of Ice Impact Loading 


The impact force on the blade arises from the momentum of the ice piece. The temporal and 
spatial distribution of this impact force is carried out with the following considerations [2J. 


a) Only the impact force component normal to the blade chord will cause the local 
deformation. 


b) Load duration, T, of impact is the time taken by the equivalent spherical ice piece (radius, 
r and relative velocity, Vr e i) to squash-up, i. e., T= 2 r/V re i- 

c) impact force is distributed using a parabolic distribution function over all finite element 
nodes within one diameter (equivalent) of the ice piece from the impact node (see 
Figures 6a and 6b). 


fe,[l-A:?/4r2)[l-y,/2r) (2) 

1 I [l-^/4r2][l- yi /2r] 

i e n 

where /; is the nodal force, f av is the component of the impact force normal to the blade 
chord, and x; and y; are the local coordinates with impact node as the origin (Figure 6a). 
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d) The impact wave travels along the chord at the same velocity as that of ice impact. Asa 
result, i th nodal loading start time, T,- , is defined by the ratio of nodal chord wise distance, 
d from the leading edge and the ice relative velocity, V re \, i.e., T,- = d/V ret , (see Figure 
6a). ' ' c 


e) Blade flexibility reduces the instantaneous ice impact loads due to reduction in the 
relative velocity. However, deflection of the blade results in ingestion of extra ice and 
therefore, there is an increase in the total momentum exchange. These factors are 
accounted and nodal forces are adjusted accordingly (see Reference [2] for details). 

Transient Response Analysis 

Transient response of the local patch is obtained by the process of modal integration. The 
modal analysis is based on the principle that statically deformed structure can be described by 
a combination of that structure's mode shapes. Once the coefficient of each mode's 
psrtidpfltion is calculated, the impact event can be simulated by the integration of a seri<?s of 
linear steps through time. The deflected shape due to ice impact is basic; therefore, only the 
first five modes of the local patch are used to compute the response. Also, damping is 
assumed to be zero as very little structural damping occurs during the impact event. 

Upon diagonalizing the governing equations of motion of undamped system [5], the 

uncoupled equations with each modal coordinate, satisfies the second order differential 
equation. 


i + $ =-j^ p i(t) ( 3 ) 

where (Op- Ki/Mj r A4,- and K{ are the i th diagonal elements of the generalized mass and 
stiffness matrices, and Pjft) is the i ^ component of the generalized force at time t. 


The general solution of the above equation, expressed in terms of arbitrary initial conditions 
a nd at t — f^,) and a convolution integral of the applied load, can be written as. 


£(»= Ht-trtti,* 


+ G(t-tn)§ /n 


+ 



G(f-r) P,(r)dt 


(4) 


where F and G, respectively are the solutions with unit displacement and unit velocity 
initial conditions. 7 

Assuming that load varies linearly between t n and t n+1 , the solution at time t n+1 can be 
obtained in closed form as (see Reference [6]), 


1322 



i,n+\ ~ F(h)£i,n + 

G(/j)£- „ + A Pi,„ + B P I>+ i 

(5a) 

>i,n+l ~~ F (W<3i,n 

+ G‘(h)Zi,n +A'Pi, n + B'Pi , n+ 1 

(5b) 


where h = t n+1 - 1 „ , F(h) = costofi; G(h) = (1/coJ situofr, F'(h) = -<*>; sincofi and G'(h) = cosafi 
for the undamped system [5]. 


For zero initial conditions (§,0 = <*,0 — 0) , the coefficients A, A , B and B are given by 


A- 1 ( «'" fek -cosmth],!)* -L-jl- 
MiO?\ °>i h I Midjfl 


sin Q)jh \ 

0)ih I 


(6a) 


A' = 


— - — sin o)ih - -Ip + 

MiiOi (Oih 


cos (Qjh 
(Oih I 


l j 1 _ cos Q)jh \ 
MiCOi\(Oih (Oih I 


(6b) 


The response at all time steps can be obtained by applying equations for &n+ 1 ^ kn+\ 

recursively. At each time step, the modal coordinates, £>, can be transformed to nodal 
displacements using modal matrix. From these displacements, large radial strains, £ of the 
leading edge elements are computed with the relation [2], 

£x= i[m^Lf (7) 

The entire calculation proceeds very quickly, due to the uncoupled nature of the modal 
responses. Among the strains at all the time steps, the maximum is picked as a 
representative quantity for the local impact damage. 


RESULTS AND DISCUSSION 


The SR-2 model unswept propfan blade [7] is chosen to study the ice impact and the 
planform of the blade is shown in Figure 7. The setting angle (orientation of the blade chord 
with respect to the plane of rotation at 75% span) is 57° and the number of blades are 8. In 
the present study, the composite material layup: Titanium skin/ Graphite-Epoxy (±45)/ 
Titanium core is used (Figure 8). For these materials, the properties are given in Table 1. 
The effective yield stress for element stiffness modification is assumed to be 206 ksi. 
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Effect of jce Size arid Ice Velocity, 


The impact in the region 50-90% of the span is known to cause severe damage to the blade 
[2]. Hence, lower and upper limits of the modeling region are taken as 50% and 90% b ade 
span for the impact analysis. For this modeling region, the impact radius is at 9.55" (75.5% 
span). These data are presented for three engine speeds, 3000, 5000 and 8000 RPM in the form 
of contour plots (Figures 9a-c). For all the engine speeds, the maximum strain corresponds 
to r = 0.8". However, the corresponding critical ice speed varies depending on the engine 
RPM. As an example, at 3000, 5000 and 8000 RPM the peaks are respectively at 100, 130 and 
190 knots. The associated peak values of the strains are 1.66%, 7.61% and 14.22%. From these 
plots, a safe range of ice impact parameter; namely, ice speed and ice size that can cause a 
specified local damage to the blade can be determined. Once the designer determines the 
tolerable strain along the leading edge of the blade, the combination of ice size and speed are 
defined by the region to the left of the level curve labelled with that strain. 

For the maximum ice size, a typical variation of strain with ice speed is shown in Figure 10. 
The strain is zero at ice speeds of 0 and 211 knots, but it reaches a maximum value at 100 

knots. This is due to the fact that only impact force normal to the chord (fsind) gives rise to 
this strain. When ice speed= 0, the impact force (/) is zero and when ice speed = 211 knots, 

the impact angle ( 6 ) is zero. Initially the value of d is equal to the stagger angle (55° at the 
impact radius) and gradually decreases to 0° as the ice speed increases (Figure 10). 

Figure 11 shows the the effect of impact radius on the strain. The strain is nearly zero when 
the impact takes place below 53% span (approx. 7.5" radius) but steeply increases as the ice 
impacts towards the blade tip. This is expected as the blade is relatively thin closer to the tip. 
Moreover, the increase in impact radius results in increase in the blade velocity and 
consequently in impact force . 

Effect of Engine Speed 

With engine RPM as the parameter (1000-10,000 RPM), the strain plot shown in Figure 10 is 
repeated in Figure 12. It can be seen that the variation of strain is smooth for engine speeds 
up to 4000 and fluctuates for the higher speeds. However, the average trend remains the 
same for all RPM and the peak strain is larger for higher engine speed. The probable cause 
of these local fluctuations is due to discrete nature of the nodal force distribution, i.e., the 
finite element nodes which are in the force distribution boundary are suddenly loaded when 
the impact wave reaches them (Figure 6a). 

Relative magnitudes of typical impact (out-of-plane) displacement distribution for throe 
engine speeds (3000, 5000 and 8000 RPM) is shown in Figure 13. The data corresponds to ice 
speed of 200 knots and ice radius of 0.8". The displacements along the leading edge remain 
the same and gradually decreases to zero toward the center of the blade. The maximum 
displacement for 5000 RPM is higher than that for 8000 RPM as the analysis ice speed (200) 
does not correspond to the critical one in both the cases. 
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CONCLUSIONS 


The ice impact on an aircraft engine blade made of layered composite is simulated. From 
the results presented in this paper, the following conclusions can be drawn. 

1. The largest ice piece results in maximum average leading edge strain. 

2. For a given engine speed, a critical ice speed exists that corresponds to the maximum 
strain. This critical speed increases with increase in the engine speed. 

3. The leading edge strain increases steeply with increase in impact radius. 

4. The strains generated from this analysis can be effectively used for structural tailoring of 
engine blades. 
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Table 1: Properties of the Blade Constituent Materials 
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Figure 2. Geometry of Ice Impact on an Engine Blade 



Figure 3. Finite Element Model for Local Impact Analysis Figure 4. Finite Element Local Coordinate System 


1327 




Figure 5. Symmetric Material Layer Configuration 
in the Composite Blade 



Figure 6a. Force Distribution Region in Local Impact Figure 6b. Parabolic Force Distribution in the Region n 
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are 10. Variation of Average Leading Edge Strain and Impact Angle with Ice Speed 
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Figure 11. Effect of Impact Location on Average Leading Edge Strain 
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Figure 12. Effect of Ice Speed on Average Leading Edge Strain with RPM as Parameter 



(a) RPM = 3000 


(b) RPM = 5000 


(c) RPM = 8000 



Figure 13. Absolute Out-of-Plane Displacements of Impact Region, Ice Size = 0.8", Ice Speed = 200 knots. 
Impact Location: 73.5% Span 
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